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ABSTRACT 

The X-ray measurements of the intra-cluster medium metalhcity are becoming more and more frequent 
due to the availabihty of powerful X-ray telescope with excellent spatial and spectral resolutions. The 
information which can be extracted from the measurements of the a-elements, like Oxygen, Magnesium 
and Silicon with respect to the Iron abundance is extremely important to better understand the stellar 
formation and its evolutionary history. In this paper we investigate possible source of bias or systematic 
effects connected to the plasma physics when recovering metal abundances from X-ray spectra. To 
do this we analyze 6 simulated galaxy clusters processed through the new version of our X-ray MAp 
Simulator (X-MAS), which allows to create mock XMM-Newton EPIC MOSl and M0S2 observations. 
By comparing the spectroscopic results inferred by the X-ray spectra to the expected values directly 
obtained from the original simulation we find that: i) the Iron is recovered with high accuracy for both 
hot (T > 3 keV) and cold (T < 2 keV) systems; at intermediate temperatures, however, we find a 
systematic overestimate which depends inversely on the number counts; ii) Oxygen is well recovered in 
cold clusters, while in hot systems the X-ray measurement may overestimate the true value by a factor 
up to 2-3; iii) Being a weak line, the measurement of Magnesium is always difficult; despite of this, for 
cold systems (i.e. with T < 2 keV) we do not find any systematic behavior, while for very hot systems 
(i.e. with T > 5 keV) the spectroscopic measurement may strongly overestimate the true value up to 
a factor of 4; iv) Silicon is well recovered for all the clusters in our sample. We investigate in detail 
the nature of the systematic effects and biases found performing XSPEC simulations. We conclude that 
they are mainly connected with the multi-temperature nature of the projected observed spectra and 
to the intrinsic limitation of the XMM-Newton-EPlC spectral resolution that does not always allow to 
disentangle among the emission lines produced by different elements. 

Subject headings: cosmology: miscellaneous - methods: numerical - galaxies: clusters: general - X-ray: 
spectrum - hydrodynamics. 



1. INTRODUCTION 

Measuring the metalhcity of the intra-cluster medium 
(ICM) represents a unique means towards a unified de- 
scription of the thermodynamical properties of the diffuse 
gas and of the past history of star formation within th e 
population of cluster galaxies (e.g., iRenzini et al.lll993fl . 
The total mass in metals is directly related to the number 
of supernovae (SNe) which exploded in the past. Fur- 
thermore, measuring the relative abundance of Iron with 
respect to a-elemcnts provides information on the relative 
fraction of type II and type la SN and, therefore, on the 
shape of the initial mass function (IMF) for star forma- 
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tion. Finally, the redshift evolution of the metal content 
of the ICM reflects the rate at which metals are released 
in the diffuse medium. This, in turn is directly related to 
the star formation history. 

From an observational point of view, measure- 
ments of the global content of different metal 
species within the ICM have been first made 



possible with the advent of the ASCA satel- 
lite fe.g.. Fukazawa et al.1 19941 iMushotzkv et al. 


199(1 iFukazawa et al. 




19981 iFinoguenov et al. 


1999; 


Finoguenov & Ponman 


119991: 'Fukazawa et al.l 


2000; 


Finogucnov et al. 200d 


Finogucnov & JonesI 


2000; 
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[Matsushita et al.ll2000HFinoguenov et al.ll200lD . Based on 
a coll ection of ASCA data for a f airly large number of clus- 
ters, iBaumgartner et al] (|2005[ ) determined the relation 
between the global ICM content of different metal species 
and the temperature. At the same time, data from the 
Beppo-SAX satellite allowed to pe rform a detailed study 
of th e spatial distribution of Iron (|De Grandi fc Molendil 
l200lh . This study showed that significant negative metal- 
licity gradients are present in "cooling flow" clusters, with 
a central meta llicity spike associated to the brightest 
cluster galaxy ()De Grandi et al.l l2004f) . More recently, 
the excellent spatial resolution of C handra has a l lowed 
to obtain Iron metallicity profiles (lEttori et afl l2002t 
Buote et al .ll2003 HDupke fc Whitdl2003l: [Humphrey et al 



20061) and maps (j Sanders et al. 1 120041 : iSanders" fc FabianI 
20061) for a significant number of clusters. Furthermore, 



the large collecting area of XMM-Newton is now allow- 
ing to study in detail the abundances of several ele- 
ments, thereb y providing a test-b e d for different models 
ofSNe yields (|de Plaa et al.ll2006l [20071: lB5hringer erall 
l2005l and reference therein), and to largely improve 
the descripti on of the metal content assoc i ated to cluster 
cool cor es (iMolendi fc Gastaldellol 12001 : iTamura et"ar 



2001 



[Ba 

[Gas t 
iMat si 



Bohringer et al 



2002 



2002 : iGas taldcUo fc Molendil 20021: 



20021: [Mat sushita et all 120031: 'Bohringer e t all 



Finoguenov et al 



Sakcllio u et al. 
~l200l 



Tamura et al...2004; IWerner et al.i ,2006.: Matsushita et al.l 
2007aD . Chandra and XMM-Newton allowed to a nalyze 



in det ail the Iron bias already noticed with ASCA (jBuotd 
l2000aD . This bias is arising when fitting with a single tem- 
perature model a plasma characterized by a temperature 
gradient or by a multi-temperat ure structure with the 
colder component belo w 1 keV (jMolendi fc Gastaldellol 
120011: iBuote et al.ll2003[ ). More recently, the Suzaku tele- 
scope is providing accurate detections of Oxygen and 
Magnesium lines and has already provided measurements 
of the distribution of g- elements in the For nax cluster 
(iMatsushita et all l2007"bh and in Abell 1060 (|Sato et al.l 
[200717 out to fairly large cluster-centric distances. Finally, 
archival studies based on Chandra and XMM-Newton ob- 
servations are now providing information on the evolution 
of the global Iron metal content out to z ~ 1.3 and show a 
decrease of about 50 per cent of the ICM metallicity from 
the present time to that redshift (jBalestra et al.l [20071 : 
iMaughan etalll2007l ). 

From the theoretical point of view, a number of meth- 
ods have been developed to describe the IGM enrichment 
associated with the cosmologica l hiera r chical build up of 



galaxy clusters. De Lucia et al.l (|2004[ ). [Nagashima et al.[ 
(|2005t ) and lBertone et al.l ([20071 ) studied the ICM enrich- 



ment by coupling TV— body simulations to semi-analytical 
models (SAMs) of galaxy formation. This approach, al- 
though quite flexible to explore the space of parameters 
relevant for star formation and galaxy evolution, does 
not provide detailed information on the spatial distribu- 
tion of metals within clusters and on the role of gas- 
dynamical processes (e.g., ram-pressure stripping, tur- 
bulent diffusion, etc.) in distr ibuting met als. In order 
to overcome these limitations, [Coral ( 20061) developed a 
hybrid technique, which is based on applying SAMs to 
non-radiative hydrodynamical simulations. Another ap- 
proach is that of simulating clusters with hydrodynami- 



cal codes, which includ ed a treatment of spe cific effects, 
such as galacti c winds ([Schindler et al.[ [20051 ) . ram pres- 
sure st ripping ("Domai nko et al.|[2006( ) and active galactic 
nuclei ([Moll et al . 200^ although coupled with simpli- 
fied descriptions of star formation. Finally, cosmological 
hydrodynamical simulations of galaxy clusters have been 
presented in the last few years, in which the production 
of metals is self-consistently related to the process of gas 
cooling and star formation, by accounting in detail for 
the contribution from dif ferent stellar populations and the 
corresponding li fe-times (iLia et al."2002';'Valdarnini 20031: 
Tornatore et al.|[20 04: Romeo ct al. 2006; Tornato re et al.l 



20071) . 



These numerical descriptions of the ICM enrichment 
have reached a good enough accuracy, in terms of physical 
processes included and resolution, so as to allow a detailed 
comparison with observational data to be performed. It 
is however clear that for this comparison to be meaning- 
ful, one has to be sure to understand observational biases, 
related to the instrumental response of the X-ray detec- 
tors onboard of the different satellites. A typical example 
of such biases in the study of the ICM is represented by 
the m easure of the X-ray temperature. [Mazzotta et al.[ 
(|2004j ) showed that the combination of a complex thermal 
structure of the ICM and of the finite energy band, within 
which X-ray spectra are measured, can bias low the value 
of the spectroscopic temperature (see also , Vikhlinin.2006[ ). 
Therefore, it is worth asking whether similar biases also af- 
fect the determination of the ICM metallicity through the 
fitting of the X-ray spectrum with single-temperature and 
single-metallicity plasma models. 

In this paper we address the issue of the biases in the 
X-ray measurement of the ICM metallicity by perform- 
ing mock observations of hydrodynamical simulations of 
galaxy clusters, which include a detailed description of the 
process of chemical enrichme nt and follow the produ ction 
of different chemical species (|Tornatore et al.[|2007D . To 
this purpose, we generate XMM-Newton observations of 
both EPIC-MOSl and EPIC-M0S2 cameras, performed 
with a new version of the X-ray MAp Simulator (X-MAS, 
[Gardini et al.|[20dih package, which now also includes the 
contribution from metal lines in the generation of the X- 
ray spectra. Taking advantage of this tool, we can care- 
fully check for the presence of systematics in the X-ray 
measurement of the metallicity, also by distinguishing the 
behavior of Iron and of a-elements, when observations are 
performed at the spectral resolution of XMM-Newton. In 
particularly we investigated the biases related only to the 
physics of the plasma. The biases related to the uncertain- 
ties on the calibration and/or on the background and/or 
on plasma model are out of the scope of this paper. As 
such, this paper mainly focusses on quantifying the pos- 
sible presence of these systematics and understand their 
relation with the complex thermal and chemical structure 
of the ICM. We defer to a forthcoming analysis a com- 
parison between results of hydrodynamical simulations of 
ICM enrichment and observational data (Fabjan et al. in 
preparation). 

The plan of the paper is as follows. The new version of 
our X-ray MAp Simulator, X-MAS2, will be presented in 
Sect.[2j The characteristics of the hydrodynamical simula- 
tions are described in Sect. [3| The procedure followed for 
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the X-ray analysis is presented in Sect.[31 Our results are 
described in Sect. [3 where we also provide a critical discus- 
sion about the effect of observing X-ray spectra with the 
typical XMM-Newton energy resolution (Sect. [5T6|) . These 
results are discussed in Sect.[6l while the main conclusions 
of our analysis are summarized in Sect. [71 

All errors quoted are at 1 ct level (68.3 per cent level of 
confidence for one interesting parameter). 

2. X-MAS2 

The numerical code X-MAS (iGardini et all l2004 f) pro- 
duces mock Chandra and XMM-Newton observations 
starting from outputs of hydro-N-body simulations of 
galaxy clusters. The software is composed by two main 
parts: the first section computes the emissivity of each 
gas particle contained in the selected field of view and 
projects it into the sky plane; the second one convolutes 
the flux reaching the telescope using the response of a 
given detector. One of the main characteristics of X-MAS 
is that it is able to create the events file containing all the 
spectral information for each "incoming simulated pho- 
ton". Moreover it properly reproduces the spatial fea- 
tures and the characteristics of the Chandra and XMM- 
Newton detectors. Furthermore the header of the final 
event file stores all the fundamental keywords common to 
true observed files. In this way the application of the rou- 
tines and analysis programs commonly used by observers 
is straightforward. In previous work the code has been 
used only in the Chandra ACIS-S3 mode to address several 
issues, like the difference between the emi ssion- weighted 
and spectroscopic temperature definitions (jGardini et al.l 
|2004| ). t he validity of the spec troscopic-like formula pro- 
posed bv lMazzotta et al] (|2004D and the possib le sources of 
bias i n the X-ray estimates of cluster masses asia et al.l 
120051 I2006D . To do that the I CM metallicity has been 
kept fixed at 0.3 solar f assuming I Anders fc Grevesselll989[ ) 
and the emissivity has been compu ted simpl y adopting 
a ME KAL model ()Mewe et al.lll985l . 1986: Li edahl et all 
1199 51. For the purpose of this paper a much more de- 
tailed analysis of the metal content of the ICM is required. 
Therefore, X-MAS has been subject to a substantial up- 
grading which affected both the first and second parts. 

2.1. First part: implementing the metallicity 

We refer to IGardini et al.l (|2004[ ) for a detailed deriva- 
tion of the incoming photons fiux, FJ . Here we report its 
final expression as a function of redshift z, luminosity dis- 
tance di. Emission Measure EM = Jy UenndV' (where 
the integrals of the electron. Tie, and hydrogen, uh, densi- 
ties are computed over the volume V) and power coefficient 

Notice that the power coefficient is a function of both tem- 
perature and metallicity. To take this latter into account, 
the metal content is split in different elements and their 
contribution is computed separately. 

Thus, the total power coefficient per each particle is 
given by: 

^;(i+.)(7^.^tot) = PV,)(r,i/)+ (2) 



where the sum is over each chemical element, Z^, con- 
sidered in the simulation (He included). The power co- 
efficients are computed assuming a single-temperature 
VMEKAL model, where we fix to 1 solar value the con- 
tribution of each element, namely He, C, N, O, Mg, Si 
and Fe. The spectrum computed by VMEKAL for each 
element, P^i^i^z)^"^' ^i), will be the sum of the continuum, 
given by the power coefficient corresponding to hydrogen, 
-Pj(i+z)(^j plus the contribution of the different el- 
ements. The factor m{Zi) represents the weight equal 
to the Zi-element mass in units of i ts Solar values (de- 
fined using the abundances reported in lAnders fc Grevessd 
|1989() . Hereafter we will compute the emissivity by us- 
ing a MEKAL model, nevertheless the code includes a fiag 
to switch to the APEC model (and cor responding VAPEC 
for multi- metal elements treatment, see ISmith et al.ll200lD 
whenever desired. 

In the new version of the simulator we have 
also modified the routine b ased on the SPH kernel 
(jMonaghan fc Lattanziol Il985l ) to distribute the particle 
emissivity, and assign the interesting quantities to the pro- 
jected pixels. The change consists in computing the inte- 
gration over the area through the products of two integrals 
along the two directions of the plane, where each integral 
has been analytically calculated at priori. This trick sig- 
nificantly speeds up the computational time. 

As in the previous version of X-MAS, we also include the 
effect of the Galactic HI absorption. Chosen a value of the 
column density Nh, we multiply the flux in each energy 
channel by an ab sorption coefficient computed adopting 
the WABS model ([Morrison k McCammonl[T983h . 

2.2. Second part: simulating EPIC observations 

Once the projected fiux has been obtained by the first 
unit of X-MAS, the second unit takes care of the simula- 
tion of cluster "observations" . This second part has been 
implemented to simulate XMM-Newton observations made 
with the EPIC camera. This is done by generating a set 
of photons, given the incoming ICM and chosen an expo- 
sure time. These are then applied to a ray-tracing pro- 
cedure designed to mimic the main characteristics of the 
telescope optical pa ths and detector responses. (See also 
iBourdin et aI]|2004D . 

The incoming set of photons {p{k, I, e)}, having position 
on mirrors [fc, I] and individual energy e, is generated from 
random realizations of the total energy expected assuming 
a given ICM flux FJ, a telescope area A, and an exposure 
time t. At this aim, we project the energies t x A x FJ 
along the line of sight, and we store the obtained quan- 
tity in a three-dimensional array (where two dimensions 
represent the position on the mirrors and the third one is 
for energy). Assuming a Poisson statistics, we can finally 
associate a discrete number of photons to the resulting to- 
tal energies. Initially we consider an array having angular 
and spectral resolutions of Ar = 3 arcsec and Ae = 20 
eV,, which are below the expected performances of the 
XMM-Newton-EPlC instruments. 

The ray-tracing procedure emulating the EPIC observa- 
tions comprises a set of filtering and redistribution func- 
tions of the photon energies and positions across the de- 
tector planes. Information about the different instrumen- 
tal effects are coming from ground and on-ffight calibra- 
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tions, and essentially provided in the XMM-Newton Cur- 
rent Calibration Files (CCF) and EPIC Redistribution 
Matrix Files (RMF). The main steps of our procedure are 
as follows: 

• redistribution of the photon coordinates, [k,l], ac- 
cording to the instrument Point Spread Functions 
(PSF), given as a function of the position on mirrors 
and energy of the incidental photons; 

• rejection of part of the incoming photons due to the 
spatially dependent cross-section of mirrors, the so- 
called "vignetting"; 

• suppression of photons falling outside CCDs or 
within dead pixels; 

• rejection of part of the incoming photons due to en- 
ergy dependent effective area of mirrors, filter ab- 
sorption, and quantum efficiency of detectors; 

• redistribution of photon energies e, according to the 
detector response; 

• addition of particles induced false detections, with 
random drawing performed following the relative 
energy spectrum expected for a given exposure 
time. 

We finally obtain a list of events, {ev{k,l,e)}, corre- 
sponding to photon detections - but also including some 
false ones due to particles -, with registered positions on 
the detector plane, [fc,Z], and energies, e. 

3. THE SIMULATED CLUSTERS 

The six simulated clusters analyzed in this paper have 
been selected from a lar ger ensemb l e of s imulations of 
19 clusters presented by ISaro et all (|2006l ): we refer to 
that paper for more details. These c lusters are extracted 
from a parent DM-only simulation (lYoshida et al.|[200l[ ) 
with a box size of 479/i"^Mpc of a flat ACDM model 
with r^om = 0.3 for the present matter density parame- 
ter, h — 0.7 for the Hubble constant in units of 100 km 
s~^Mpc~^, (78 = 0.9 for the r.m.s. fluctuation within a 
top-hat sphere of 8 h~^Mpc radius and ilob = 0.04 for the 
baryon density parameter. Mass resolution is increased 
inside the interesting regions by using th e Zoomed Ini- 
tial C ondition (ZIC) technique proposed bv lTormen et al.l 
(|l997t ). Besides the low-frequency modes, which were 
taken from the initial conditions of the parent simula- 
tion, the contribution of the newly sampled high-frequency 
modes was also added. Once initial conditions are created, 
we split particles in the high-resolution region into a DM 
and a gas component, whose mass ratio is set to reproduce 
the assumed cosmic baryon fraction. In the high-resolution 
region, the masses of the DM and gas particles are set to 
muM = 1-13 X IO^/i-IMq and mgas = 1.7 x IO^/i-^Mq, 
respectively. The Plummer-equivalent softening length 
for the gravitational force is set to epi — 5ft-~^kpc, kept 
fixed in physical units from z = 5 to z = 0, while being 
epi = 30ft.^^kpc in comoving units at higher redshifts. 

The simulations have been ca rried out with a ver- 
sion of the GADGET-2 code (|SDringe]| I2005D . which 

13 http://www.MPA-Garchmg.MPG.DE/gadget/ 



includes a detailed trea tment of chemical enr i chmen t 
from stellar evolution (jTornatore et all |2004 |2007[ ). 
GADGET-2 is a parallel Tree+SPH code with fully 
adaptive time-stepping, which includes an integration 
scheme which explicitly co nserves energy and entropy 
(|Springel fc HernquistI 12002 ), the effect of a unifor m and 
evolving UV background ( Haardt &: Madaul Il996[ ) , star 
formation from a multiphase interstellar medium and a 
pres cription for galactic winds t riggered by SN explosions 
fsee lSpringel fc Hernauistjr2003l for a detailed description), 
and a numerical scheme to suppress artificial viscosity far 
from the shock regions (see Dolag et al. 2005). In the orig- 
inal version of the code, energy feedback and global metal- 
licity were produced only by SNII under the instantaneous- 
recycling approximation (IRA) . We have suitably modified 
the simulation code, so as to correctly account for the life- 
times of different stellar populations, to follow metal pro- 
duction from both SNIa, SNII, as well as from low and 
intermediate mass stars, while self-consistently introduc- 
ing the dependence of the c ooling function on metallicity 
(jSutherland fc Dopitalll993l ) . A detailed description of the 
im plementation of t hese algorithms is presented in detail 
by Tornatore et al.l (|2007T ) . The simulations analyzed here 



assume the power-la w shape for initi al stellar mass func- 
tion, as proposed by ISalpeteil (|l955t) and galactic ejects 
with a speed of 500 km s~^. 

The main global properties of the simulated clusters are 
reported in Table [TJ In the sample there are two big ob- 
jects with mass (computed inside the radius i?2oo which 
comprises an over-density of 200 over the Universe criti- 
cal density) larger than IO^^Mq, a medium-mass cluster 
(Af2oo — 3.51O"'^'*M0) and three smaller systems, namely 
G676, G914 and G1542. Since these last three objects 
have very similar characteristics and behaviors, we chose 
to report our results only for the first one, as example. 
In Fig. [T] and Fig. [2] we show the projected maps for the 
spectroscopic-like temperature and for the metallicity, re- 
spectively. The map size is equal to the field of view of 
XMM-Newton, i.e. 30 arcmin: since we put our clusters 
at redshift z = 0.06, this corresponds to 2088 kpc, using 
the cosmological model here assumed. The different metal 
maps refer to Iron and Oxygen (representative of the a- 
elements). From these maps it appears that the largest 
clusters, Gl.a and G51, have a rather complex structure, 
while the smaller cluster G676 appears as more relaxed, 
with Gl.b representing an intermediate case. 

In the following we briefly describe the dynamical state 
of our simulated clusters. This is fundamental to explain 
the sample characteristics which could influence the pres- 
ence of any bias. The two massive clusters, Gl.a and G51, 
both have an active dynamical state and strong inhomo- 
geneities in the temperature maps (upper panels of Fig. [T]) . 
They show a monotonically increasing mass history, espe- 
cially lately: this means that there is a significant amount 
of material which continuously in-falls inside the cluster 
virial region. All these small structures can be recognized 
both in the temperature maps (Fig. [T|) and in the metal- 
licity maps (Fig. [J) of these objects (see also the photon 
images reported in Fig. [3]). Those are cold blobs which 
have not had enough time to thermalize in the hot ICM. 
Their presence does not affect the results of our analy- 



X-MAS2: Metallicity Systematics 



5 



sis, since they have been identified and excluded from the 
X-ray analysis (the masked regions are shown as green cir- 
cles in Fig. [Tl see Sect. S]). The metal maps show rather 
elongated structures, which trace the direction of the most 
recent merger events. 

The cluster Gl.b is close to Gl.a, lying in the same re- 
simulated Lagrangian region. It has experienced two con- 
secutive important major merger events and it has not had 
enough time to relax. This clearly appears from the com- 
plex features in the temperature map (see also the presence 
on the photon image of a bright extended source on the 
West side of the centre). The core of the merged halo is 
detectable by eye from the metal maps. Quite interest- 
ingly, the Oxygen distribution can be used to reconstruct 
the merging dynamics. 

The cluster G676 is an isolated system, which presented 
a relevant merging event at z « 0.7, but after then it has 
had time to relax. The object can be considered in dy- 
namical equilibrium at present time, in the sense that its 
global velocity dispersion and its mass are almost constant 
in time after the merging event. Its homogeneous temper- 
ature map, round flux contours and spherical metallicity 
distribution confirm this picture. Similar considerations 
can be done for the other two small objects, G914 and 
G1542 (not shown in the figures). 

4. THE X-RAY ANALYSIS 

The simulated clusters have been processed through X- 
MAS2 in order to produce mock observations with the 
EPIC-MOSl and EPIC-M0S2 cameras and with an ex- 
posure time of 200 ksec. As mentioned before, the field of 
view is 30 arcmin, corresponding to 2088 kpc in physical 
units. The images have been created by integrating along 
the line of sight a box of 10 Mpc centered on the clusters. 
In our analysis we assume that XMM-Newton is perfectly 
calibrated. 

4.1. Spatial analysis 

We create the photon images in the [0.7-2.0] keV band 
and use them to detect the presence of compact X-ray 
bright gas regions. These corresponds to cool and very 
dense gas clumps present in the simulation and likely pro- 
duced by overcooling effects. These clumps are generally 
not observed in real clusters so to avoid contamination we 
remove these regions from our analysis. For this purpose 
we used the wavelet de composition algorithm proposed by 
IVikhlinin et al.l (|l998l ). Parameters are settled to local- 
ize sources with a peak above 4.5tT and to identify the 
surrounding region above the 2.5a level. The selected ar- 
eas are shown as green circles in the temperature maps 
(Fig.[T|). In Fig. [3] we show the photons images, extracted 
in the [0.7-2] keV energy band, corrected by vignetting, 
out-of-time events, point spread function, background and 
out-field-of-view component, binned to 3.2 arcsec. 

4.2. Spectral analysis 

The spectral analysis has been performed by extract- 
ing the spectra from the event file in concentric annuli, 
centered on the brightest X-ray peak of the cluster. The 
first annulus has an internal radius of 50 kpc (the cen- 
tral region is excluded from the analysis) and an external 



radius which is 1.25 times larger. All the remaining an- 
nuli are logarithmically equi-spaced with the same step of 
1.25. During the spectrum extraction, we also subtracted 
all photons coming from the masked regions, already men- 
tioned in the previous section, and from box regions con- 
taining the CCD gaps (clearly detectable as horizontal and 
vertical lines in the photon images displayed in Fig. [J]). 
The spectra have been fitted using th e C statistics in the 
XSPEC 11.3.2 package (|Arnaudl [l996l ) and following two 
procedures widely used by observers: 

- Method 1. We fit the data with an absorbed 
VMEKAL model in the [0.4-8] keV energy band, 
leaving at the same time as free all the parameters 
of interest: temperature, Iron, Silicon, Magnesium, 
Oxygen and normalization. This approach is faster 
and more automatic but can lead to artificially bi- 
ased results because of the degeneracy of the pa- 
rameters. 

- Method 2. This is a four steps procedure, (i) we 
fit the data with an absorbed MEKAL model in 
the [0.4-8] keV band to obtain the temperature; 
metallicity and normalization are considered free 
parameters; (ii) we fix the temperature and use a 
VMEKAL model in the same energy band to re- 
cover the Iron abundance; the other three metals 
are left free; (iii) we keep frozen temperature and 
Iron to measure the Oxygen in the [0.4-1.5] keV 
band and (iv) finally we freeze the values of tem- 
perature, Iron, Oxygen to estimate the Magnesium 
and Silicon abundances in the [1.2 3.2] keV band. 

In both fitting methods we fixed the Galactic absorption 
Nh and the redshift to the input values used to simulate 
the observations with X-MAS2 (i.e. Nh = 5 x 10^° cm'^ 
and z = 0.06). 

5. RESULTS 

From the spectral analysis we recovered the projected 
profiles for temperature. Iron, Silicon, Magnesium and 
Oxygen and we compared them to the values directly ex- 
tracted from the hydrodynamical simulations. The sim- 
ulation profiles are computed projecting the quantities 
along the line of sight, within the same radial bins used 
in X-ray analysis and masking out the same regions de- 
scribed in Sect.m The temperatur e is computed using th e 
spectroscopic-like definition, Tsl (jMazzotta et al.l [20041 ) . 
summing over all particles with temperature higher than 
0.5 keV. For the metal profiles we use instead the emission- 
weighted formula: 

Z,„=J2Z^XW,/J2W^, (3) 

where Wi is the emissivity of the i-th particle computed in 
the [0.4-8] keV energy band considering the contribution 
of each particle metals. 

The results for Gl.a, G51, Gl.b, and G676 are pre- 
sented in Fig. m The black lines refer to the simulated 
input values, while spectroscopic values determined from 
the full-band-fitting (Method 1) and the ones recovered in 
the narrow band are shown by the green and red lines, 
respectively. The vertical error bars are at la level while 
horizontal bars represent the width of the radial bin. 
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5.1. Temperature 

The recovered spectroscopic temperature matches weh 
the input temperature leading to a maximum relative er- 
ror of 5 per cent for hot systems. This difference is ex- 
pected for objects with a complex temperature structure 
(see Fig.[T]) and steep temperature gradients ( see Fig.|4]), as 
already delineated in iMazzotta et all (12004 ). It is worth 
noting that the excellent agreement is also valid for the 
three cooler systems, for which the temperature falls from 
2 keV at 2 arcmin 150 kpc) down to 1 keV at 6 arcmin 
{'^ 450 kpc). Indeed, this result is expected since these 
objects show a regular and homogenous temperature map 
and any weighting (i.e., mass- weighted, emission- weighted 
and spectroscopic-like) give the sam e consistent resu lt in 
the limit of isothermal systems [Mazzotta et al.|[2004l) . 

As expected, we do not find any difference between the 
spectroscopic temperatures recovered using Methods 1 and 
2, i.e. between MEKAL and VMEKAL fitting. In fact we 
find that for the most massive clusters the two estimates 
almost perfectly coincide; in the other cases the differences 
are always smaller than 2 per cent. 

5.2. Iron 

Thanks to the powerful emission of both Fe-L and Fe-K 
lines, the spectroscopic determination of the Iron abun- 
dance is recovered with great accuracy and presents rel- 
atively small error bars. The strength of Fe-K lines for 
objects having a temperature larger than 3 keV and the 
one of Fe-L lines for systems with temperature lower 
than 2 keV (see Fig. [5]) ensures an appreciable agree- 
ment for both cold and hot clusters. In each radial bin 
the Iron abundance derived spectroscopically matches the 
emission- weighted value from the original hydrodynamical 
simulation at 3tT level. Moreover the relative difference, 
[Fespcc-Fesim]/Fesim, is less than 5 per cent. The only ex- 
ception is for Gl.b, the 2-3 keV cluster (see Fig. 2]). For 
this object we detect a systematic overestimate with re- 
spect to the input profile which can be as high as 20 per 
cent. The reason of this discrepancy - related to the com- 
plex temperature structure of this cluster (already evident 
in the corresponding panel of Fig. [1]) and to its particular 
temperature range - will be extensively discussed later in 
Sect. 16.21 Some spatial bins of G676 show an underesti- 
mate of Iron, which is between 10% and 30%. This fea - 
ture is caused by the Iron bias described in lBuotd ((20003) • 
Among our cold clusters, G676 is the only one presenting 
this characteristic; we will further investigate this aspect 
in Sect.O 

It is worth noting that the Iron spectroscopic abun- 
dances obtained from Method 1 and Method 2 are not 
significantly different. In fact, the two procedures differ 
just by keeping fixed the temperature in the second case, 
but the possible degeneracy between Iron and Tempera- 
ture is weak because of the large number of counts we had 
(larger than lO''). The maximum discrepancy between the 
two methods is smaller than 10 per cent and it is present 
only for some radial bins of the coolest objects. 

5.3. Oxygen 

The spectroscopic measurements of Oxygen, as well as 
Magnesium and Silicon, are more uncertain in the inner 
regions of the clusters. The effect is larger for the hottest 



systems, Gl.a and G51, and decreases for objects with 
smaller masses. The explanation resides on projection ef- 
fects: in fact the spectrum of the central region is ob- 
tained by integrating along the whole line-of-sight, thus it 
includes contributions from plasma at temperatures which 
can significantly differ in the case of most massive clusters, 
because of their steeper radial profile. For the systems 
with temperature lower than 3 keV, the spectroscopic de- 
termination of the Oxygen abundance always agrees with 
the input values inside l-2a. On the other hand, for the 
largest clusters, Gl.a and G51, the Oxygen detection is 
systematically overestimated, even if the large error bars 
reduce the statistical significance in most radial bins. In- 
creasing the exposure time, i.e. the photon number, re- 
duces the error and the systematic effect becomes more 
evident. To confirm this result, in Fig. [5] we report the 
spectroscopic result for Gl.a, our hottest cluster, for which 
we simulated a 1 Msec exposure. From this plot we notice 
how the spectroscopic values significantly overestimate the 
input value by a factor of 3. 

As a general behavior, we notice that Oxygen obtained 
by fitting the spectra with Methods 1 and 2 are similar for 
the smallest objects and with discrepancies always inside 
the error bars. Nevertheless, there is a tendency of obtain- 
ing an abundances larger than the input one when Oxy- 
gen is measured using the broad band (Method 1 shown 
as green dashed points in Fig. 2]), being the difference in- 
creasing with the temperature. 

5.4. Magnesium 

Generally speaking, the measurement of Magnesium is 
very difficult because, at the spectral resolution of XMM- 
Newton, its line energy is well within the energy range of 
the Fe-L line complex (see Sect. 15.61 and Fig. In few 
cases it is really impossible to detect any Magnesium con- 
tent, and when a detection is obtained we get always very 
large errors. Things improves slightly if we use Method 2 
but not by much. 

Magnesium lines are stronger for cool systems: in fact, 
in the cases of Gl.b and G676 (see the third and fourth 
columns of Fig. [3]), we are able to trace the overall pro- 
file, even if no strong constraint can be obtained. In any 
case for these cold clusters we do not find any systematic 
bias. On the contrary, the 6 keV cluster, G51, shows large 
discrepancies in the external bins where the spectroscopic 
estimates are between 200 and 400 per cent higher than the 
input values (see Fig. 0]). Finally, for Gl.a no conclusion 
can be reached since its high temperature (all radial bins 
have values above 8 keV) makes the Magnesium extremely 
difficult to be detected (see also Fig. [5]). 

5.5. Silicon 

The difference between the results obtained using Meth- 
ods 1 and 2 is significant also in the measurements of the 
Silicon abundances. Shrinking the energy band improves 
the estimates, mainly for temperature larger than 2-3 keV. 
For the cluster Gl.b, the choice of a narrow energy band 
(red lines in Fig. [4]) helps in preventing the underestimate 
of this metal by 40-50 per cent. G51 and Gl.a show a sim- 
ilar tendency, even if this result is less significant because 
of the larger error bars. The colder clusters, instead, do 
not show any particular difference between the two spec- 
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tral methods: in both cases the input abundance of Sihcon 
is well recovered from the X-ray analysis at 1-1.5 a level. 

In general the agreement between the X-ray results and 
the input values from the hydrodynamical simulations is 
remarkably good at all temperatures. It is important to 
say, however, that the silicon line energy is located in re- 
gion where, due to a combination of CCD and mirrors, the 
effective area falls down rapidly and after 2 keV it shows 
a strong edge due to Au and Ir. For real observations this 
may be a problem as this spectral region of the detector is 
very difficult to calibrate. However, this problem will not 
be present for high-redshift clusters, for which the Silicon 
line is redshifted towards the softer region. 

5.6. Thermal Breemstrahlung spectra as seen by 
XMM-Newton 

Our main goal is to investigate how well we can measure 
the metal abundances with the availability of the XMM- 
Newton EPIC camera spectra. In order to better under- 
stand the results of our analysis we illustrate here some 
general features of thermal breemstrahlung spectra ob- 
served with the same spectral resolution of XMM-Newton 
detectors. We generate with XSPEC six spectra at different 
values of temperature: 1, 2, 3, 5, 8, and 10 keV. These are 
convolved with the instrumental response of the MOS-1 
camera (almost identical results can be obtained consider- 
ing MOS-2). The adopted plasma model is VMEKAL, for 
which we set to zero the abundances of all metals, except 
O, Mg, Si and Fe, which are fixed to their solar values. The 
resulting spectra are shown as red lines in FigO In the 
same plot different color lines refer to the separate contri- 
bution of each metal: Oxygen (green) , Magnesium (blue) , 
Silicon (cyan) and Iron (black). Decreasing the plasma 
temperature, we can notice that the emission lines become 
more evident and the ratio between Fe-L (Iron lines around 

1 keV), and Fe-K (Iron lines around 6 keV), rapidly in- 
creases. The Oxygen lines, due to both OVII and OVIII 
lines, show two bumps around 0.6-0.8 keV which are more 
evident for systems with temperature lower than 3 keV. 
The Oxygen lines are located very close to the Fe-L com- 
plex, whose tail in the soft energy region, given the spectral 
resolution of XMM-Newton, can interfere with Oxygen de- 
termination for systems with temperature around 1-2 keV. 
At temperatures larger than 3 keV, the Oxygen and Fe- 
L features are well separated. The Magnesium lines are 
quite strong only for very cold systems with temperature 
around 1-3 keV. Nevertheless, also in this case, the spec- 
tral resolution of the XMM-Newton EPIC-MOS cameras 
is not good enough to allow to distinguish this line from 
the Fe-L complex, which has an extremely powerful emis- 
sion in this temperature range. The importance of the 
Magnesium lines is rapidly decreasing. Silicon, instead is 
present in a more isolated region of the spectrum (close to 

2 keV), therefore this guarantees a good measurement of 
its abundance for systems with temperature up to approx- 
imately 8 keV. Another important property of the energy 
region around 2 keV is that different temperature spec- 
tra with the same normalization have there the same flux. 
This means that the Silicon measurement is not signifi- 
cantly influenced by the presence of a multi-temperature 
plasma or, in other words, by the uncertainties related to 
the continuum determination. 



It is important to stress that the possibility of accurately 
determining the emission lines is subject not only to the 
mutual interaction of different elements, but also to the ca- 
pability of measuring with accuracy the Iron lines (which 
can strongly infiuence Oxygen and Magnesium) and the 
continuum. 

6. DISCUSSION 

6.1. Differences between the methods used in the spectral 

analysis 

Considering the results presented in Sect. [5] (see also 
Fig- SI), we can conclude that the two strategies followed 
in the X-ray analysis give almost identical estimates in the 
case of temperature and Iron. The only relevant difference 
is that in Method 2 we freeze the temperature to compute 
the Iron. The effect is to eliminate one free parameter, re- 
ducing the degeneracy, which, however, is extremely weak 
in our case, because of the large number of counts per 
each radial bin. The most evident differences between the 
two procedures are for the abundance measurements of the 
a-metals (Oxygen, Silicon and Magnesium): in all cases 
we obtain a better result when we fix the continuum and 
the Iron contribution and we fit the spectrum in a nar- 
row band centered on the corresponding set of lines. We 
made several tests to individuate the optimal energy band 
to adopt. At the assumed redshift, i.e. z = 0.06, the best 
choices are [0.4-1.5] keV for Oxygen and [1.2-3.2] keV for 
both Magnesium and Silicon. Notice that the bands need 
to be large enough to include part of the continuum to bet- 
ter estimate the line emission and to evaluate the effects 
of each single element on the continuum. In particular 
this shrewdness is very important for Oxygen, for which 
the continuum contribution is relevant on a wide energy 
range (see the green line in the first panel of Fig. [5]). Using 
Method 2 instead of Method 1 allows to avoid an overes- 
timated result for Oxygen. Moreover it permits to detect 
Magnesium which would not be recovered without freezing 
the Iron. Finally it increases the spectroscopic determina- 
tion of Silicon avoiding a possible underestimate. 

6.2. Iron overestimate for systems at 2-3 keV 

The Iron content of the ICM is perfectly recovered for 
hot systems, Gl.a and G51, because their Fe-K lines are 
extremely well determined (see Fig.|3]and the last 3 panels 
of Fig. [S]). Iron is also accurately estimated for cold sys- 
tems (only exception are some bins of G676 in Fig. [4]- see 
Sect. 16. 3p . since their Fe-L lines are extremely strong (see 
the first two panels of Fig. [5]). The only system for which 
we obtain a systematic overestimate in all radial bins is 
Gl.b, the 2-3 keV object. Notice that this discrepancy 
is significant, as evident from the error bars in the corre- 
sponding panel of Fig. |4l The problem is due to the fact 
that the Gl.b temperature is very close to the one where 
there is a transition between the relative importance of 
the lines (Fe-L or Fe-K) used in the determination of the 
global Iron content. Moreover, its temperature structure 
is quite complex: consequently inside the same radial bin 
we are averaging different temperatures which sometimes 
are characterized by a large Fe-K contribution, while other 
times show a strong Fe-L. Summing over spectra with both 
high Fe-L and high Fe-K has the net effect of increasing 
the final amount of Iron. This effect can be explained by 
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looking at Fig. \7\ In the first panel we combine together a 
plasma at 1 keV (green line) with a plasma at 4 keV (dark 
line): the resulting spectrum is shown by the red line. It 
is evident that the first plasma contributes to bump the 
Fe-L lines, while the second one has the role of bursting 
the Fe-K lines in the resulting global spectrum, which is 
then characterized by an Iron content which is larger than 
what expected considering its temperature. 

In order to quantify this effect we perform a test with 
a likely composition of the gas (see the lower panel of the 
same figure). The input data is the spectrum resulting 
from the combination of two plasmas: the first one with 
temperature Ti = 2 keV, metallicity Zi = 0.2 and normal- 
ization Ki = 1; the second one with T2 = 3 keV, Z2 = 0.1 
and K2 = I- We produce few spectra with difi^erent num- 
ber counts and fit them with a C statistics. The results 
are reported in Table [3] They show a perfect consistency 
for the temperature, but a significant overestimate of the 
Iron spectral value with respect to the emission-weighted 
one. The discrepancy is highly significant and depends on 
the number counts: for 4.500 counts the disagreement is of 
almost 40%. Increasing the statistic number the difference 
is alleviated but not solved. Even observing the cluster for 
2 X 10^ sec and thus having 9 x 10^ counts (which is 1- 
2 orders of magnitude above the typical observed counts 
rate, see e.g. iBalestra et aLl l2007^. the discrepancy is still 
present and of order of 15%. The spectrum corresponding 
to the best-fit (red line) of the 4.5 x 10"'^ counts spectrum 
is compared to the input data in the same panel of Fig|71 
while in the bottom box we show the corresponding resid- 
uals. Notice that in the residual behavior there are not 
evident features which can make us suspicious that the fit 
is not good or that it can be affected by some bias on the 
Iron lines. 

The implications of these results are extremely impor- 
tant. Iron is in fact the metal to which all the other a- 
elements are referring, thus overestimating its content can 
lead to wrong conclusions on the past stellar population. 

This effect may also explain, at least partially, the re- 
cent results obtained bv lBaumgartner et al.l (|2005f ). After 
stacking in temperature bins 273 clusters observed with 
ASCA , they studied the global metal abundances as a 
function of the mean temperature. They found that in all 
temperature bins (from 0.5 keV to 12 keV) the emission- 
weighted Iron abundance within the cluster region selected 
to contain as much flux as possible is around 0.2-0.25 (in 
solar units), with the only exception of the two bins be- 
tween 2 and 4 keV, where a global Iron value of 0.43-0.47 
Zq is measured. 

6.3. Fe bias for cold systems 

Another bias affecting the estimate of the Iron con- 
tent is present in some bins of the cold system G676 (see 
Fig. 3]). The system atic underesti mate is know in litera- 
ture as "Iron bias" (jBuotd [200031 10) . However, this bias 
gives at maximum an underestimate of order of 20% and 
it is present only in G676, but it does not influence the 
other two cold systems. The poor evidence of this bias in 
our sample can be explained as follow. 

The Iron bias introduced by iBuotd ()2000aD is caused 
by the fact of pretending to fit with a single temperature 
model a plasma which is, instead, characterized either by a 



combination of different temperatures (multi-temperature 
plasma) or by a strong temperature gradient. Observa- 
tionally, the bias is seen to affect the central cooling core 
regions where typically there is a positive temperature 
gradient {dT/dr > 0) as well as both negative metallic- 
ity gradient [dZ/dr < 0) and negative emission gradient 
{dEM/dr < 0). The combination of these factors leads 
to enlarge the bump of Fe-L shell (the colder plasma with 
higher metallicity and emission excites the lines on the 
soft part of Fe-L shell). In our work, we are not in this 
regime. On the contrary, both temperature and metallicity 
gradient are shallow inside the spatial bins we considered. 
Moreover the temperature gradient in the our external re- 
gions is going in the opposite direction (dT/dr < 0) with 
respect its behavior in the cooling regions. 

Another reason for which we do not detect a strong Fe 
bias in our system is because it decreases if the spectral 
fitting is made in an energy band sufficiently large to de- 
termine thBron^muu^ or, in other wor ds, the temp e rature 
profile (|Buotell2000bl ) with accuracy. iBuote et al.l (|2003D 
suggested to extend the minimum of the fitting energy 
band down below 0.6 keV in order to limit the effect of 
this bias. With our analysis we respect this demand since 
the minimum energy considered is 0.4 keV. In addition, a 
robust determination of the continuum is guaranteed by 
the high signal to noise that we have up to 4 keV. The 
fact that a good Iron measurement depends on a correct 
temperature measurement is clear comparing in Fig. 2] the 
G676 temperature and Iron panels. 

6.4. Oxygen and Magnesium for hot clusters 

In Sect. [5] we showed that the abundance estimates for 
both Oxygen and Magnesium do not show any problem of 
systematics for systems colder than 2 keV. On the other 
hand, for the hot clusters, like Gl.a and G51, we found 
a clear bias producing a significant overestimate. Notice 
that these elements are the best a-element indicators, thus 
they are a fundamental tool for studying the stellar evolu- 
tion and formation history through the ratio between SNe 
la and SNe II. It is therefore necessary to understand in 
detail the origin of the discrepancy. At this aim we inves- 
tigate a number of possible reasons which depend on the 
ICM physics implemented in the original hydrodynamic 
simulation, on the complexity of the plasma temperature 
and on the features of the X-ray spectra. 

• Presence of small undetected cold blobs. 

Clusters Gl.a and G51 show in their photon im- 
ages (Fig. [3]) a large number of compact gas clouds 
spread in the ICM which are usually associated 
to cold blobs in the temperature maps (Fig. [J). 
Most of them have been detected by the wavelet 
decomposition algorithm and then excluded from 
our analysis, but some of them are still present, in 
particular in the external regions. These structures, 
which have a very low temperature, could in prin- 
ciple be responsible for an overestimate for Oxygen 
or Magnesium, since their importance is greater for 
cold systems. To investigate this problem in more 
detail we generate new event files for Gl.a and G51, 
but using only the particles having a temperature 
larger of 5 keV: in this way we are sure to avoid any 
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contamination from cold (diffuse and clumped) gas. 
To these event files we apply the whole procedure 
describe above. As expected, using the wavelength 
algorithm, we do not find any compact cool source 
and we just exclude the second brightest blob, lo- 
cated on the East side of G51 centre. The com- 
parison between the spectroscopic results and the 
simulation profiles of Oxygen for Gl.a and Magne- 
sium for G51 (re-computed excluding the particles 
with temperature smaller than 5 keV) are shown 
in Fig. [51 it is clear that the discrepancies are not 
alleviated. 

• Dynamical state of the cluster. Since the 
bias is evident only for the two largest clusters 
which clearly have a complex morphological struc- 
ture, while it does not affect the coolest systems 
which are more symmetric and regular, its relevance 
could be enhanced by the presence of thermal in- 
homogeinities in the ICM. However, we notice that 
Gl.b does not show this phenomenon, even if it rep- 
resents a perturbed object. To further investigate 
this possibility, we generate another series of events 
files fixing the temperature of all the particles at 
a given constant value (8 keV) but letting all the 
metals equal to their original values. In this way 
we create an object with a completely flat tempera- 
ture structure. We repeat the analysis on these new 
images, finding again the same discrepancy for the 
Oxygen and Magnesium abundances. 

• Fe-L over Fe-K ratio. The spectra we have 
extracted from the images are the sum of con- 
tributions coming from the plasma present along 
the line of sight and having different temperatures 
and metallicities. The resulting spectra, therefore, 
could have a Fe-L over Fe-K ratio different from the 
typical value reported by the VMEKAL model in 
correspondence of the measured temperature. The 
Iron content of hot systems is estimated basically 
using the Fe-K systems which are strongly more 
powerful than the Fe-L lines at those temperatures. 
The Fe-L group of the spectra can be, thus, weaker 
than what we expect from a VMEKAL model where 
we estimate the Iron content by fitting the Fe-K 
lines. In this situation one way to "fill the gap" 
between the true Fe-L of the spectra and the ex- 
pected one of the model is to increase the emission 
power by other elements present in the same energy 
band. We performed a detailed comparison of the 
spectral Fe-L and Fe-K measurements and the cor- 
responding model expectations. We conclude that 
the Magnesium estimate can be affected for this rea- 
son. In fact for G51, which is a system showing a 
large disagreement for Mg lines, the Iron determi- 
nation from the Fe-L is always higher than what we 
derive from the Fe-K lines. On the other hand, this 
behavior is not present in Gl.a and thus cannot ex- 
plain the disagreement found for Oxygen, which is 
also displaced by Fe-L at temperatures larger than 
8 keV (see Fig.E]). 

• Continuum determination. 



Oxygen is a weak line at temperatures larger than 5 
keV, thus a wrong determination of the continuum 
can strongly influence the measurement of this el- 
ement. In particular an underestimate of the con- 
tinuum can originate an overestimate of the Oxy- 
gen abundance. To have a rough estimate of this 
effect we simulate a spectrum as sum of a 6 keV 
plasma with Oxygen equal to 0.1 (in solar units) 
and a 8 keV plasma with Oxygen fixed to 0.2. The 
two plasmas have been weighted in the same way 
by assigning a unity normalization to both. The 
fake spectrum has been convolved with the XMM- 
Newton response and simulated assuming an expo- 
sure time of 1 Msec. This large integration time al- 
lows us to have few millions of counts per spectra. 
We performed a first fit using a statistics and 
a VMEKAL model where we let free temperature. 
Oxygen and normalization: the resulting temper- 
ature agrees with the spectroscopic-like estimate. 
This temperature has been then fixed in the fol- 
lowing tests, where we also fixed the normalization, 
but to different values, starting the input theoreti- 
cal value of 2 and then decreasing it slightly every 
time. In this way we are able to study the effect that 
only a wrong determination of the continuum has 
on the spectroscopic determination of Oxygen. The 
expected emission-weighted abundance for Oxygen 
is 0.149. This value is perfectly recovered when we 
force the normalization to be the theoretical one. 
A result compatible at la level is still obtained us- 
ing a lower normalization, up to 1.995. Decreasing 
further the normalization by 0.5 or 0.75 per cent, 
the spectroscopic determination of Oxygen becomes 
much larger than the expected one. This simple 
test shows us how Oxygen is very dependent on the 
goodness of the continuum determination and sug- 
gests that the measurements done through XMM- 
Newton in high-temperature systems can be signif- 
icantly biased. 

7. CONCLUSIONS 

The Intracluster Medium represents a perfect laboratory 
to study the physics of all the cluster subcomponents, like 
its galaxies and the stars in galaxies. In particular the ICM 
metallicity is a rich source of information to test the stellar 
evolution models. In this work we deeply investigated all 
possible systematic biases present in the metallicity mea- 
surements from X-ray spectra having the same spectral 
resolution of XMM-Newton. At this aim we improved the 
X-ray MAp Simulator in order to correctly consider the 
metal information which is provided by our simulations. 
The six clusters here analyzed cover a wide range of tem- 
perature and have different dynamical state. This allows 
us to determine how the biases we analyzed are effectively 
influenced by the complexity of the thermal structure and 
by the value of the temperature itself. 

Performing the X-ray analysis we tested two standard 
procedures adopted by X-ray observers: first, we fitted the 
spectra in a large band leaving free at the same time all 
the parameters of interest (T, Fe, O, Mg, Si); second, we 
recover the temperature and the Iron content from a fit 
in a large band and then we measure the a-elements, O, 
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Mg, and Si, in a narrower band centered on the lines. Our 
main results are shown in Fig. |4] and can be summarized 
as follows: 

- Using narrow bands to measure all the elements 
beside Iron is in general better for all clusters, in- 
dependently of their temperature. This procedure 
prevents to overestimate Oxygen and to underesti- 
mate Silicon. Furthermore it increases the proba- 
bility to detect the Magnesium content. 

- Iron is perfectly recovered as soon as either Fe-K 
or Fe-L are the dominating lines: this occurs when 
the plasma has a temperature larger than 3 keV or 
smaller than 2 keV, respectively. In our sample we 
have an object showing a complex map temperature 
(Fig. [1]) with values ranging between 2-3 keV (see 
its temperature profile in Fig. 3]) . Simulating a sim- 
ple model which describes this system in XSPEC we 
found that the spectroscopic Iron is overestimated 
by up to 15-40 per cent depending on the num- 
ber statistic: decreasing the source counts leads to 
an increase of the disagreement. This result has 
strong implications. In fact one of the main goals 
of measuring the ICM metallicity is the derivation 
of the [a/Fe] and [Si/Fe] ratios (see the Introduc- 
tion). Therefore, the fact that we detect a clear 
systematic effect in the Iron measurements can af- 
fect the results that we can derive from them about 
the stellar formation history. 

- Oxygen is well measured for clusters with temper- 
ature lower than 3 keV, while we find an evidence 
of overestimate for the hottest system in our sam- 
ple, which has a temperature profile ranging from 8 
to 12 keV (see Fig. We check carefully some of 
the possible reasons that can originate this effect. 
The explanation came out to be the problematic 
balancing between the determination of the contin- 
uum and the weakness of this line at temperature 
higher than 8 keV. A discrepancy smaller than 1 per 
cent on the normalization is sufficient to explain the 
Oxygen overestimate. 

- The measurement of Magnesium is difficult at all 
temperatures: for cold objects (with T < 2 keV) its 
line is, in fact, lying in the same spectral range of 
strong Fe-L lines and thus its detection is hidden, 
while at high temperatures it is a very weak line. 
We found a systematic overestimate of Magnesium 
in hot clusters: this effect can be related to the fact 
that the [Fc-L/Fc-K] ratio is higher in the spectra 
than in the model used to fit them. Thus to fill 
the Fe-L group lines there is the necessity to pump 
the emission of other elements present in the same 
energy band, that in our case, is Magnesium. 

- Finally, from the X-ray analysis of the spectra we 
can measure with good precision Silicon for all the 
objects in our sample independently of their dynam- 
ical state. This is due to the fact that the value of 

http:/ /projects. iasf-roma.inaf.it/edge/ 
http:/ /constellation. gsfc.nasa.gov/index. html 



the spectra around 2 keV is only slightly dependent 
on the temperature. 

The method used in this paper to investigate the metal- 
licity distribution takes an enormous advantage by the 
use of high-resolution simulations where the hydrodynamic 
physics is treated including many different physical pro- 
cesses which are relevant for the thermal and chemical 
evolution of the ICM. Other aspects which can play a sig- 
nificant role, like AGNs, cosmic rays and magnetic field, 
are not included in these simulations but could be added 
in the next future. At the present the state of our sim- 
ulations does not fully reproduce the cluster core, where 
more of the X-ray information comes from, but it is a 
prope r description of th e temperature gradient outside 0.1 
i?200 (jPratt et al.ll2007[ ). Nevertheless, it is worth stress- 
ing that our results and the corresponding explanations 
are not dependent on the kind of simulations used. In 
fact, when we found some disagreement between the input 
profiles from our simulation and the X-ray spectroscopic 
measurements, we investigated the problem by analyzing 
the spectral properties of fake plasma spectra generated 
using XSPEC . For this reason our results can be con- 
sidered valid for observations of spectra generated by a 
plasma with a multi-component temperature. The only 
restrictions of our results are due to the fact that we fo- 
cused our effort on a fixed spectral energy resolution, like 
the XMM-Newton one. The powerful of the analysis of the 
ICM metal content can radically change in the next fu- 
ture with an X-ray satellite having a micro-calorimeter on 
board. EDGE (Explorer of Diffuse Emission and Gamma- 
ray burst Explosions) Ell and Constellation- are two 
projects that have been proposed with the intent of in- 
corporating such technology, that will allow to solve the 
blending of Fe-L and Oxygen lines. 
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Fig. 1. — Spectroscopic-like temperature maps for the mock clusters: Gl.a (top-right panel), G51 (top-left panel), Gl.b (bottom-right 
panel), G676 (bottom-left panel). The color scale is linear, ranging from a minimum value of 0.1 keV up to a maximum value, which is equal 
to 13, 11, 4.5, 3 keV for Gl.a, G51, Gl.b and G676, respectively. We overplot (black curves) the flux contours (level equi-spaced by a factor 
of 2) and the compact X-ray bright cool regions detected using the wavelet decomposition algorithm (green circles). 



Table 1 

The main physical properties of the simulated clusters analyzed in this paper: R20O AND iisoo IN UNITS OF Mpc/h 
(the RADIUS OF THE SPHERES WITH A DENSITY 200 AND 500, RESPECTIVELY, TIMES THE CRITICAL DENSITY OF THE 

Universe); the masses inside these radii, M200 and M500, in units of and the mass-weighted 

TEMPERATURES IN KEV, rMW,200- 





^200 


-R500 


M200 


-^^500 


?MW,200 


Gl.a 


1.771 


1.202 


12.92 


10.09 


7.37 


G51 


1.675 


1.120 


10.92 


8.16 


6.44 


Gl.b 


1.151 


0.740 


3.55 


2.36 


2.63 


G676 


0.727 


0.490 


0.89 


0.68 


1.42 


G914 


0.718 


0.485 


0.86 


0.66 


1.39 


G1542 


0.725 


0.484 


0.88 


0.66 


1.32 
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Fig. 3. — Photons images in the [0.7-2] keV energy band of the clusters in fig [T] Gl.a (top-right panel), G51 (top-left panel), Gl.b 
(bottom-right panel), G676 (bottom-left penal). The images are binned to 3.2 arcsec, background subtracted, and corrected for vignetting 
and out-of-time events 
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Fig. 4. — Comparison of the spectroscopic and projected input radial profiles. From top to bottom we show the Temperature, Iron, Oxygen, 
Magnesium, Silicon profiles, respectively. Each column correspond to a different mock cluster (from left right to right Gl.a, G51, Gl.b and 
G676 clusters). Solid black line represents the spectroscopic-like temperature and the emission- weighted metal abundance. Green dashed and 
solid red points with relative error bars refer to the spectroscopic values derived by using Methods 1 and 2, respectively (see text). 
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Fig. 5. — The red lines represent the spectra obtained from XSPEC reshifted at 2 = assuming a VMEKAL model with all the metals 
equal to zero but O, Si, Mg and Fe which are fixed to the solar values: the temperature is varying from 1, 2, 3 (top panels starting from the 
left) to 5, 8, 10 keV (bottom panels). The separated contributions of the different elements are also plotted: Oxygen, Magnesium, Silicon 
and Iron are shown by green, blue, cyan and black to lines. All spectra are convolved with the response of EPIC MOS-1 camera; very similar 
results can be obtained considering the EPIC MOS-2 camera. 
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Fig. 6. — The radial profile of Oxygen in Gl.a observed with an exposure time of 1 Ms. The meaning of the lines is as in Fig. l4l 
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Fig. 7. — Top panel: Spectra of a plasma at Ti=l keV (green line), T2 =4 keV (black line) and their combination (red line). Bottom panel: 
a fake spectrum with 450.000 counts obtained by combining a plasma with Ti = 2 keV and Z\ = 0.2 and a plasma with T2 = 3 keV and 
Zi = 0.1. The red line shows the resulting best-fitting model (see text for more details). The lowest box presents the corresponding residuals. 



Table 2 

Values of Oxygen abundances (with their la error bars), derived by fitting a spectra, which is the combination 
OF two plasmas at 6 AND 8 keV (more details in the text). In the fitting procedure we fix the temperature 
at the spectroscopic-like value and then we fix also the normalization, considering the its theoretical value 

(2) and smaller values (1.995, 1.990, 1.985). The percentages (within brackets) represent the difference 
between the values recovered from the spectral analysis and the expected emission- weighted one: 0=0.149. 



2.000 1.995 1.990 1.985 

Oinin 0.126 0.148 0.175 (15%) 0.192(29%) 

O 0.141 0.162 (9%) 0.184 (23%) 0.198(38%) 

Omax 0.155 0.176 (18%) 0.198 (33%) 0.220(48%) 

x\ed (70 d.o.f.) 1.08 1.08 1.12 1.20 
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Fig. 8. — The radial profiles for Oxygen in Gl.a (top panel) and Magnesium in G51 (bottom panel), obtained considering only particles 
with temperature larger than 5 keV. Symbols and lines are as in Fig. |4] 



Table 3 

Values of Temperature and Iron abundances (with their la error bars), derived by fitting a spectra of a T=2 

keV and Z=0.2 plasma combined with another plasma of T=3 keV and Z=0.1. The percentages (within 
brackets) represent the difference between the values recovered from the spectral analysis and the expected 
emission-weighted one: Fe=0.144. The expected spectroscopic-like value of the temperature is equal to 2.45. 
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